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Artificial quasicrystals are nowadays routinely manufactured, yet only two naturally-occurring 
examples are known. We present a class of systems with the potential to be realised both artificially 
and in nature, in which the lowest energy state is a one-dimensional quasicrystal. These systems are 
based on incommensurately charge-ordered materials, in which the quasicrystalline phase competes 
with the formation of a regular array of discommensurations as a way of interpolating between 
incommensurate charge order at high temperatures and commensurate order at low temperatures. 
The nonlocal correlations characteristic of the quasicrystalline state emerge from a free energy 
contribution localized in reciprocal space. We present a theoretical phase diagram showing that the 
required material properties for the appearance of such a ground state allow for one-dimensional 
qnasicrystals to form in real materials. The result is a potentially wide class of one-dimensional 
qnasicrystals. 

PACS numbers: 71.45.Lr,61.44.Br,61.44.Fw 


While many quasicrystals have been artificially created 
since their first synthesis in 1984, only two naturally- 
occurring quasicrystals have been identified to date - 
both found in the same Siberian meteorite [TH3]- In 
the quasiperiodic, never-repeating, pattern of two or 
more unit cells defining a quasicrystal (QC) [H |5], lo¬ 
cally changing a single cell may require adjustment of an 
infinite number of other cells in order to maintain the 
quasiperiodicity [5]. This global property makes it non¬ 
trivial to grow a QC cell-by-cell. Proposed mechanisms 
in two or higher dimensions include recognising vertex 
matching rules which imply ‘forced tiles’, random tiling 
models, and relaxation processes from non-quasiperiodic 
systems [7]. 

In (quasi) one-dimensional (ID) systems, QCs consist 
of a quasiperiodic sequence of long and short unit cells. 
The global nature of ID QCs is apparent from their cor¬ 
respondence to projections of regular 2D crystals 00. 
One way to experimentally access properties of ID QCs is 
to consider instead a setup containing two periodic sub¬ 
systems that are incommensurate with respect to one an¬ 
other. The sequence of subsystem species encountered as 
one traverses the system is aperiodic, and can be mapped 
onto a similar aperiodic or quasicrystalline sequence of 
unit cells (see Fig. [^a)). Such incommensurate systems 
have recently been employed in a number of experimen¬ 
tal studies of QCs: optical waveguide arrays simulating 
ID crystals with incommensurate periodic on-site poten¬ 
tials; cold atom condensates in optical lattices; Moire 
superlattices in graphene derivatives; and in theoretical 
descriptions of incommensurate charge density waves 0- 
na. While these are of much interest in their own right, 
they only constitute quasicrystals by way of mathemat¬ 
ical analogy. One key difference with true quasicrys¬ 
tals for example, is that incommensurate systems have 


no lower bound on the spatial separation between sites 
of the two subsystems; it is only the sequence which is 
quasiperiodic (see also the Supplementary Material). 

In this Letter we describe a mechanism for generating 
ID quasicrystals based on coupling two incommensurate 
subsystems and allowing them to adjust to minimize their 
overall free energy. We take as a specific example the case 
of incommensurate charge density waves, showing that, 
under specific conditions, it is energetically favourable 
for the atomic lattice to adopt a quasiperiodic structure. 
The resulting state is thus a true quasicrystal, transcend¬ 
ing the mere mathematical equivalence between ID QCs 
and incommensurate charge density modulations on top 
of a periodic atomic background. We show that the re¬ 
quired material properties are physically realistic. Fi¬ 
nally, we outline the criteria needed for this mechanism 
to apply more generally. 

Low-dimensional materials are prone to develop charge 
order. In the ideal ID case, a Peierls instability always 
yields a spontaneous periodic modulation of the electron 
density, accompanied by periodic displacements in the 
atomic lattice with the same wavevector |13j . In real 
quasi-lD (or higher-dimensional) materials, charge order 
typically occurs only if the electron-phonon coupling is 
sufficiently strong [HI [15]. In that case the order arises 
when the electronic susceptibility overcomes the cost of 
populating phonon modes, so that: 

/dk|5(k,k-fq)|2/^fi2i:2^^>ScCq (1) 

J ^k+q — Ck 

where (/(k, k -1- q) is the electron-phonon coupling, which 
is generically momentum- and orbital-dependent, /(ck) 
is the Fermi-Dirac distribution for an electron with mo¬ 
mentum k and energy Ck, and HtOq is the energy of the 
phonon mode which softens in the charge ordering tran- 
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sition. The full susceptibility x(q), defined by the left- 
hand side of Eq. Q, has its maximum at q, which deter¬ 
mines the wavevector Q(Tcdw) = q at which the charge 
order and atomic displacements first form. In the pres¬ 
ence of strong nesting, guaranteed in ID, q equals the 
nesting vector. In higher dimensions it is determined by 
the momentum dependences of both the electron-phonon 
coupling and the electronic dispersion [15]. 

The evolution of the charge density’s wavevector Q(T), 
as temperature is lowered beyond the onset temperature 
Tcdwj can be modeled by a Ginzburg-Landau expansion 
of the free energy naiizj: 



The modulation on top of the average electron density is 
written as p(x) = ■i/i cos((/)(x)), where for a density wave 
without defects (^(x) = Q • x. For notational convenience 
the quadratic part of the free energy has been written in 
reciprocal space, where it is proportional to the full elec¬ 
tronic susceptibility y. This term favors the formation of 
charge order at wavevector q. The other terms are more 
conveniently written in real space. Those proportional 
to c„, with n > 2, represent the local coupling between 
the atomic lattice and the charge modulations. These 
favor charge modulations locally commensurate with the 
lattice, by giving a non-zero contribution to F whenever 
n(p{K) = K • X, with K a reciprocal lattice vector of the 
atomic lattice mM- 

At the transition temperature TcdWi the order param¬ 
eter amplitude tp is small, so that the lower-order terms 
in Eq. ([^ dominate, and the order forms at wavevec¬ 
tor Q(Tcdw) = q- As temperature is lowered, ip in¬ 
creases, and the terms proportional to c„ begin to com¬ 
pete with the second-order term. Since the coupling to 
the lattice has an effect only for the parts of the density 
wave that are locally commensurate, the charge order 
with ^(x) = Q • X will not gain energy from the final 
terms for any incommensurate value of Q. Instead, the 
charge-ordered state can lower its energy by adopting 
a locally commensurate structure (/>(x) = K • x — i5(x), 
where (5(x) introduces a regular array of broadened phase 
slips (discommensurations) [I9j . which render the aver¬ 
age propagation vector incommensurate even if the local 
structure is predominantly commensurate |18j . The evo¬ 
lution as temperature decreases further is characterized 
by a sharpening of the initially-broad discommensura¬ 
tions [mds]. 

In many materials, Q(T) starts out incommensurate 
at the onset of the charge-ordered state, evolves towards 
a close-by commensurate value, and locks in at a second 
transition temperature Tic below which the charge or¬ 
der remains commensurate. Well-known examples can be 
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FIG. 1. (color online) Possible charge orderings when the 
electron density (blue) develops an instability with wavevector 
Q incommensurate with the ion lattice (pink). Here we used 
a wavelength oi 2 tt/Q = y/ll.7 ~ 3.42 lattice spacings. (a) 
Perfectly incommensurate charge order. Labeling the blue 
peaks A and the pink ions B, the sequence of A’s and B’s is 
a quasiperiodic tiling of the units ABBB and ABBBB m- 

(b) Discommensuration state, which is locally commensurate 
with period 3, but maintains the average wavevector Q PE). 

(c) A quasicrystal composed of a quasiperiodic arrangement 
of period 3 and period 4 unit cells, maintaining the average 
wavevector Q. This state is reached by shifting each peak A 
onto the nearest ion B. Dotted lines show the incommensurate 
state for comparison. In all cases the ion positions will adjust 
to the altered electronic charge distribution (not shown). An 
alternative comparison between these states is provided in the 
Supplementary Material. 


seen in 2iL-TaSe2 and TTF-TNQ In other mate¬ 

rials, the evolution of Q(T') starts out in a similar man¬ 
ner, but remains incommensurate to zero temperature. 
Examples include Kq.sMoOs, TbTea, and NbSe 2 [2TH23] . 
A proliferation of discommensurations can describe ei¬ 
ther type of behavior, depending on the value of the ratio 
c„/x(q) in the free energy of Eq. Q, which determines 
the strength of the coupling to the lattice relative to the 
maximum strength of the full susceptibility. The pres¬ 
ence of discommensurations in these materials has been 
experimentally verified, for example in scanning tunnel¬ 
ing microscopy experiments on NbSe 2 [24] . 

The discommensuration state takes advantage of all 
the terms in the free energy of Eq. ([^. Nevertheless, one 
can imagine conditions under which it may not be the 
optimal configuration of (/>(x). Consider incommensu¬ 
rate charge order which forms at Tcdw with an incom¬ 
mensurate wavelength A = 27r/g lying close to halfway 
between three and four lattice spacings (Fig. Ba)). In 
that case, it is not a priori obvious whether the coupling 
to the lattice will prefer the charge order to lock in at 
period three or four. Choosing one, and introducing dis¬ 
commensurations such that the local electronic density 
modulations are commensurate while the average prop¬ 
agation vector remains incommensurate, necessarily re¬ 
sults in a high density of rapid phase changes (shown in 
Fig-Bb)). These sharp discommensurations involve high 
momentum components of p(q) (higher harmonics of q), 
and are therefore energetically costly with respect to the 
quadratic term in F. An alternative way of distorting 
the incommensurate charge order and gaining local lock- 
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FIG. 2. (color online) Calculated X-ray diffraction pattern for 
the discommensuration state (blue) and quasicrystal (pink) 
shown in Fig. (purple indicates overlapping peaks). The 
black line depicts the approximation to the full susceptibility 
X used to calculate the phase diagram of Fig. Both types 
of order have their main peak at q, but the quasicrystal has a 
dense covering of secondary peaks close to the main diffraction 
spot, which benefit from the nonzero x in that region. The 
units are arbitrary, but both diffraction patterns have the 
same scale. 

in energy, while retaining the overall incommensurability, 
is depicted in Fig. [^c). This configuration is obtained 
by moving every maximum in the charge density wave 
directly onto the atomic position closest to it. The re¬ 
sult is a charge modulation that is locally commensurate 
everywhere, but with two local periodicities. The locally- 
commensurate patches are smoothly connected (although 
higher derivatives are not smooth), and therefore do not 
involve the high-momentum components characteristic of 
the discommensurations. For a more detailed description 
of the difference between this state and the discommen¬ 
suration state, see the Supplementary Material. 

It is clear from the incommensurability that the pat¬ 
tern of period three and four modulations in the multiply- 
commensurate structure never repeats itself. In fact it 
can be shown, using a cut-and-project construction, that 
the arrangement is quasiperiodic |15j . In reality the ions 
will be displaced towards the charge maxima simultane¬ 
ously with the displacement of the charge maxima to¬ 
wards the ions. The structure shown in Fig. 0^) will 
therefore not just contain incommensurate charge modu¬ 
lations or a quasiperiodic electronic structure, but will 
constitute a true ID quasicrystal. At the onset tem¬ 
perature of the multiply-commensurate order, there will 
be a crystal-to-quasicrystal phase transition breaking all 
translational symmetries of the original lattice. 

The X-ray diffraction pattern of a quasicrystal is char¬ 
acterized by a dense array of sharp peaks 0 [7j. As 
shown in Fig. the Fourier transform of the multiply- 
commensurate charge order indeed contains such an ar¬ 
rangement, centered around the average incommensurate 
wavevector q. Real quasi-ID materials have susceptibil¬ 
ities x(‘l) which decay over a finite range of momenta. 
The clustered peaks around q can therefore also gain en¬ 
ergy from the quadratic term in F. This situation should 
be contrasted with that of the Fourier transform of the 



FIG. 3. (color online) (a) the predicted phase diagram emerg¬ 
ing from Eq. for Q = 27r/vTL7 and fe = 1. (b) the free 
energies of the different ordered states as a function of 'ipca/a 
for fixed ‘tpdlcs = 1.2/8. The abbreviations used are: QC- 
quasicrystal, DCn-discommensuration state with commensu¬ 
rate areas of period n, ICDW-incommensurate charge density 
wave, and Q = 1/n-commensurate charge order with period 
n. 

discommensuration state, also shown in Fig. In that 
case, although the central peak gains energy, the higher 
harmonics typically lie far away in fc-space, and cannot 
utilize the non-zero width of the susceptibility. The dif¬ 
ference in the extent to which these structures profit from 
the quadratic term in F can render the quasicrystal en¬ 
ergetically favourable to the discommensuration state. 

To quantify the relative stabilities of the quasicrystal 
and discommensuration states, we consider again the free 
energy of Eq. <§■ We model the full susceptibility as 
x(q) = a/((T^(|q| — q)^ + 1), where the width 1/a of the 
peak is taken to be small but non-zero. The ratio b/a 
of the quartic and quadratic component of F at q = q 
determines the overall magnitude of the order parame¬ 
ter ijj, but does not affect the form of 4>{'x.). Whether 
the discommensuration or quasicrystalline state is favor¬ 
able is thus determined by the ratios Ci)/cs) and 
(V’^C 3 )/(V’^a). In Fig. 0a) we present the phase dia¬ 
gram as a function of these two ratios. The energy of the 
quasicrystalline state is compared to that of an array of 
discommensurations with optimized density and widths. 
A comparison of the energies of individual states, as a 
function of ipcz/a for fixed ipc^/cs, is shown in Fig. 0b). 

The energy of the incommensurate state does not de¬ 
pend on the coupling to the lattice, and is constant in 
ipcs/a. The commensurate state on the other hand does 
not gain any energy from the susceptibility term, but 
takes full advantage of the lock-in term, and decreases 
linearly in energy as ipc^/a increases. At high coupling 
to the lattice the commensurate state thus becomes fa¬ 
vorable to the incommensurate. The discommensuration 
state interpolates between these two extremes, having its 
main Fourier component at the incommensurate value, 
but being locally commensurate. It takes advantage of 
all terms in the free energy, and outperforms both the 
commensurate and incommensurate states. 

The quasicrystal also has its dominant Fourier com¬ 
ponent at the incommensurate value. At intermediate 










4 


values of ipcs/a it additionally takes advantage of the 
spread in x(q) and only pays a modest price for not be¬ 
ing optimally commensurate everywhere. In this regime, 
the quasicrystal is favored over the discommensuration 
state, giving rise to the extended area of stability shown 
in Fig. [^a). At either extreme of '^ca/a, however, the 
discommensurations win out: the spectral weight of the 
dominant QC Fourier peak is constant, so does not ap¬ 
proach the weight of the incommensurate state in the 
limit of zero lattice coupling. The discommensuration 
state is therefore favorable to the quasicrystal at very 
low values of tpc^/a. For strong coupling to the lattice, 
the locally-commensurate nature of the quasicrystal al¬ 
lows it to gain lock-in energy, with the overall energy gain 
a weighted average of the C 3 and C 4 terms. However, if 
these two coefficients are not equal, the quasicrystal can¬ 
not gain as much energy as the discommensuration state, 
which is always locally commensurate with the most fa¬ 
vorable commensurability. The discommensuration state 
therefore also wins out over the quasicrystal at very high 
values of the coupling to the lattice. 

Based on these considerations, the prerequisites for 
this type of ID quasicrystal to emerge in any real (nat¬ 
ural or artificial) material are: (i) a quasi-lD charge- 
ordered material with an incommensurate wavevector ly¬ 
ing close to halfway between two commensurate values; 
(ii) a sharp, but not infinitely sharp, susceptibility; (hi) 
a strong, momentum-independent electron-phonon cou¬ 
pling. The latter is necessary to ensure that the coupling 
of the charge order to the lattice does not strongly favor 
any particular commensurate value of the propagation 
vector. 

Quasi-ID materials contain strongly-bonded chains of 
atoms with relatively weak inter-chain coupling. Their 
Fermi surfaces typically contain nearly-flat pieces, giv¬ 
ing a strongly-peaked susceptibility. The charge ordering 
wavevector is therefore close to 2kp, and is generically in¬ 
commensurate. The question then arises as to why we are 
not overwhelmed with naturally-occurring quasicrystals 
in quasi-lD materials. First, the charge modulation p is 
small compared to the average charge density, implying 
that a dense set of peaks neighboring q in the diffraction 
pattern may not be straightforward to distinguish exper¬ 
imentally. Second, in real materials it may be expected 
that the lock-in energy is different at different commen¬ 
surate values El, putting most materials close to the left 
side of Fig. [^a). Note, however, that the order param¬ 
eter amplitude '0 increases with decreasing temperature, 
and thus enhances the value of ipC/^/c^. Conversely, at 
the onset of ordering, ip is vanishingly small, and the 
free energy is dominated by the second-order term. This 
implies that the charge-ordered phase starts off at the ori¬ 
gin of the phase diagram, and moves upwards and to the 
right only as temperature decreases below Tcdw- Third, 
except in cases of extremely strong electronic nesting, 
the momentum- and orbital-dependence of the electron- 


phonon coupling is generically essential for charge order 
to develop in real materials El- These points aside, the 
results of Fig. suggest that a targeted search may re¬ 
veal previously-unnoticed quasicrystalline phases at low 
temperatures in quasi-ID materials with incommensu¬ 
rate charge order. Experimentally, the quasicrystalline 
charge order could be identified for example using mi¬ 
croscopy to directly identify arrangements of neighbour¬ 
ing period-four and period-three cells, or using diffraction 
probes to find a fractal distribution of sharp peaks. 

To conclude, we have shown that for an extended range 
of material parameters quasi-lD systems minimize their 
free energy by developing quasiperiodic tilings of two in¬ 
equivalent unit cells. This provides a wide class of ID 
materials with the potential to become quasicrystalline. 
Candidate materials, both natural and artificial, could 
be revealed by a directed search for materials exhibiting 
incommensurate charge order, with a wavelength close to 
halfway between two commensurate values, whose elec¬ 
tronic structures are well-nested, and whose electron- 
phonon couplings are strong but momentum indepen¬ 
dent. The quasicrystalline state also extends the set of 
possible low-energy configurations of incommensurately 
charge-ordered systems, adding to the known incommen¬ 
surate and commensurate states, and their interpolation 
via discommensurations. 

Existing proposals for quasicrystal growth mechanisms 
have centered on single-cell additions at defect sites, 
or local adjustments based on propagating defects [7]. 
The quasicrystalline state discussed here, on the other 
hand, arises from a global adjustment of two periodic 
but incommensurate parent structures. The required co¬ 
ordination for this global effect in real space derives from 
the fact that it is a localized peak of the susceptibility in 
reciprocal space which determines the dominant wavevec¬ 
tor of the quasicrystalline state. 

Since the free energy contains two incommensurate pe¬ 
riods, it is perhaps unsurprising that a ground state ex¬ 
ists which respects the lack of translational symmetry. 
Closely related states have been found in both the quan¬ 
tum Hall effect in the Tao-Thouless limit [25], and the 
formation of Wigner lattices in the Hubbard model [26| . 
In the latter case, the Wigner lattice is invoked to ex¬ 
plain the optical spectra of nearly-quasicrystalline states 
in the ID charge-ordered material TTF-TCNQ p71128] . 
Both systems, however, concern high-denominator ratio¬ 
nal commensurate periods rather than the incommen¬ 
surate situation considered here, and neither depends on 
the particular free energy arguments arising from Eq. (§. 

The mechanism for forming quasicrystals considered 
in this paper has the potential to arise whenever two 
incommensurate periodicities are simultaneously present 
in the same system. If an appropriate coupling between 
the subsystems can be introduced in the presence of a 
global constraint, a free energy similar to that of Eq. (§ 
may be expected to apply, and a quasicrystalline state 
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will generically emerge over an extensive portion of phase 
space. 
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Supplementary Material 

To bring to the fore the differences between the various 
charge-ordered phases discussed in the main text, it is 
convenient to express the modulation on top of the aver¬ 
age charge modulation in each state as: 

p(x) = ^ cos (K • X — d(x)) . (3) 

Here, K is a wave vector commensurate with the lattice, 
which we will set to K = 27r/(3a) from here on. The 
commensurate charge density wave (CCDW), incommen¬ 
surate charge density wave (ICDW), discommensuration 
state (DC), and quasicrystalline state (QC) can then be 
distinguished by considering the distinct behaviour of the 
phase d(x) in each. This is shown in Fig. below. 

For the CCDW, <5(x) is zero. For the ideal ICDW, d(x) 
is the straight line (5(x) = (Q — K)x, with Q the incom¬ 
mensurate wave vector. The discommensuration phase is 
given by sections of zero slope connected by vertical steps 
at regular intervals. The height of each step is such that 
the density wave jumps forward by precisely one lattice 
spacing. For K = 27r/(3a), the step height is thus 27r/3. 
The width of the horizontal sections is such that the av¬ 
erage slope of the entire stepped structure is equal to Q. 
The quasicrystal also has average slope Q, but consists 
of horizontal sections separated by sections with slope 
(5(x) = (K' — K)x, where K' is a second commensurate 
wave vector: in this case 27r/(4a). The sections of zero 
and non-zero slope in this phase form a quasiperiodic 
sequence. 

The DC phase often arises in real materials as the most 
energetically-favourable way of combining electronic den¬ 
sity modulations with a lattice coupling m- In prac¬ 
tice, the discommensurations in these materials will not 
be perfectly sharp. Instead, they broaden slightly and 
the connection in i5(x) with the neighbouring horizontal 
regions is made smooth. We follow McMillan in account¬ 
ing for this using a Fourier expansion with coefficients, 
controlling the discommensuration width, selected so as 
to minimize the free energy of the DC state. The QC 
phase, on the other hand, is created by starting from an 
ideal ICDW and then shifting the peaks of the electronic 
density modulation directly onto the nearest atomic po¬ 
sition. This creates a quasiperiodic sequence consisting 
of sections with one of two well-defined commensurate 
wave vectors. 

The QC phase and ICDW both have quasiperiodic el¬ 
ements. For the ICDW, the sequence of atomic positions 
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FIG. 4. (color online) Possible charge orderings discussed in 
the main text. For all cases the phase 5(x) defined in Eq. (§ 
is shown divided by tt, as a function of position in units of 
the lattice spacing a. With these conventions, the commensu¬ 
rate phases are represented by straight lines with rational val¬ 
ues for their slopes. The incommensurate phase is a straight 
line with an irrational slope. The discommensuration phase 
consists of horizontal sections connected by a regular, peri¬ 
odic array of discommensurations, which are vertical lines of 
height 2/3. The spacing of the discommensurations is such 
that the phase returns to the dashed incommensurate line at 
every discommensuration. Finally, the quasicrystal consists 
of a quasiperiodic sequence of horizontal (period 3) sections 
of width 3, and sections at the slope corresponding to com¬ 
mensurate order with period 4, with width 4. The aperiodic 
sequence can be constructed by always choosing the section 
whose end-point stays closest to the dashed incommensurate 
line. 


and charge maxima encountered as the material is tra¬ 
versed is quasiperiodic. This type of ‘quasiperiodicity’, 
however, is generic to any superposition of two incom¬ 
mensurate structures, including even things like two par¬ 
allel picket fences seen from a distance. In contrast, the 
QC phase is quasiperiodic in each of its components indi¬ 
vidually. Moreover, it also consists of precisely two unit 
cells. Each piece of period four within the QC state is 
locally indistinguishable from any other piece of period 
four. In the ICDW, the atomic neighbourhood of each 
peak in the electronic structure is perfectly unique. In 
that sense, it contains infinitely many unit cells, rather 
than two. We therefore classify the QC state as a true 
one-dimensional quasicrystal, in contrast to the ICDW, 
which is only a generic quasiperiodic sequence of charge 
maxima and atoms. 
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